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Abstract. Intergalactic Lya opacity suggests that H I was reionized at z ~ 6, 
while He II reionization was delayed to z ss 3. Both epochs are in conflict with 
inferences from CMB optical depth (WMAP) which suggest early reionization at 
z = 10 — 20. One of the major contributions of FUSE to cosmological studies has 
been the detection of He II Lya ( "Gunn-Peterson" ) absorption in the spectra 
of AGN at redshifts z > 2.03. Spectra of He II absorbers, taken in concert with 
corresponding H I (Lya) lines, allow us to fix the epoch of helium reionization at 
ZHeii ~ 2.8±0.2. Here, I review FUSE observations of He II absorption, together 
with their implications for the sources and transport of ionizing radiation over 
30-50 Mpc distances through the IGM. FUSE observations of He II absorption 
toward HE 2347-4342, combined with Keck and VLT observations of H I, are 
consistent with photoionization by QSOs, with a wide range of intrinsic spectral 
indices, and modified by filtering and reprocessing in the IGM. The He II/H I 
ratio (77) exhibits variations over 1 Mpc distance scales (Az « 10~ 3 ). Intrigu- 
ingly, this 77-ratio is also correlated with Lya filaments and voids. The ionizing 
radiation field appears to be softer (higher He II/H I) in the voids. These void 
regions may be ionized by local soft sources (dwarf starburst galaxies), or the 
QSO radiation may softened by escape from AGN cores and transport through 
denser regions in the "cosmic web". The differences in ionizing spectra may 
explain the 1.4 Gyr lag between H I and He II epochs. 



1. Overview: Reionization of the Intergalactic Medium 

Before beginning a discussion of the FUSE contributions to the study of He II 
reionization, I provide an overview of the cosmological significance of "reioniza- 
tion" . This term refers to the sudden increase in ionization fraction of hydrogen 
gas in the intergalactic medium (IGM), initiated at redshifts z > 6 by the first 
sources of ionizing radiation. After the epoch, of recombination at z rec ps 1300, 
the IGM entered a long period of "dark ages" when the IGM was mostly neutral. 
This period ended sometime before z ~ 6, perhaps as early as z = 10 — 20. 

The ionization state of the IGM affects the transmission of radiation through 
the Lya absorbers (H I and He II), as well as the electron-scattering optical 
depth to the cosmic microwave background (CMB) radiation. As individual 
ionizing sources (stars, galaxies, quasars) began emitting ionizing radiation, re- 
gions of ionized gas were produced, preceded by ionization fronts (I-fronts) in 
the IGM. These isolated ionized zones (H II regions, He III regions) eventually 
overlapped (Gnedin 2000), after which the IGM became transparent to Lya ra- 
diation. Residual absorption continued as a "picket fence" of discrete, redshifted 
Lya absorbers, often termed the "Lya forest" - see spectra in Songaila (1998). 
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The rapid increase in IGM ionization occurred as a topological phase change 
in the volume filling factor of ionized zones. The longer the first ionizing sources 
were on, the greater the volume of the zones. In order to effectively absorb all 
the Lya radiation at z = 6 requires a neutral fraction of just xhi = 2 x 10~ 4 , 
even in gas that is underdense by a factor of 10 relative to the mean. Simulations 
of this reionization process (Gnedin 2000, 2004) show that the transition from 
neutral to ionized is quite extended in time, between z = 5 — 10. The first 
stage (pre-overlap) involves the development and expansion of the first isolated 
ionizing sources. The second stage marks the overlap of the I-fronts and the 
disappearance of the last vestiges of low-density neutral gas. Finally, in the post- 
overlap stage, the remaining high-density gas is photoionized from the outside. 
The final drop in neutral fraction and increase in photon mean free path occur 
quite rapidly, over Az « 0.3 (see Fig. 1 of Gnedin 2004). 

Reionization is tied closely to the onset of galaxy formation beginning at 
z > 20. For the ionization of He II, and somewhat for H I, reionization depends 
on the rise in the numbers and luminosities of active galactic nuclei from z = 
5^2. According to numerical simulations (e.g., Ricotti, Gnedin, & Shull 
2002, herafter denoted RGS), the earliest epochs of galaxy formation began at 
z « 20 — 30 triggered by H I and H2 cooling of protogalactic halos with virial 
temperatures T v - ir = 10 3 K to 10 4 K. RGS found that the average star formation 
rate density, psfr, increased monotonically from 10 -3 to 10 _1 M & yr _1 Mpc -3 
between z = 20 and 10, according to the approximate formula: 

SFR = (0.003 M yr" 1 Mpc _3 )10 (20_2)/5 . (1) 

The effects of reionization can be seen indirectly in the electron-scattering 
optical depth of the cosmic microwave background (CMB), as we discuss in 
the next section. However, the most direct observations of the IGM ionization 
state come from their influence on the Lya line-blanketing opacity of the IGM at 
wavelengths shortward of the H I and He II Lya lines at (1215.670 A)(l + z) and 
(303.782 A)(l + z), respectively. As the amounts of H I or He II increase, these 
line-blanketed regions grow into black ( "Gunn- Peterson" ) absorption troughs, 
whose presence can be used to identify high-redshift galaxies (Steidel et al. 
2002). The best probes of the epoch of reionization have come from analyzing 
the weak transmission of AGN continuum flux in these troughs, as the high- 
redshift IGM transforms from a neutral to ionized medium. The He II absorption 
becomes detectable in the far ultraviolet at redshifts z > 2.03 with FUSE and 
at z > 2.79 with HST. From the latest Lya data, the reionization epochs occur 
at z ~ 6.1 i 0.3 for H I (Becker et al. 2001; Fan et al. 2002; Gnedin 2004) and 
at z » 2.8 ±0.2 for He II (Kriss et al. 2001; Smette et al. 2002; Shull et al. 2004; 
Zheng et al. 2004). 

2. Electron-Scattering Optical Depth of the IGM 

An indirect measure of reionization can be extracted from an analysis of fluctua- 
tions in the CMB radiation, particularly the large-angle temperature-polarization 
correlations. By combining the first-year data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) with other results on the CMB, Lya forest, and 
galaxy large-scale structure, Spergel et al. (2003) found a high value of the 
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CMB optical depth to electron scattering, r e = 0.17 ± 0.04. This optical depth 
implies a surprisingly large redshift of early IGM reionization, ranging from 
11 < z r < 30 at 95% confidence (Kogut et al. 2003). More recent estimates of r e 
and other cosmological parameters (Tegmark et al. 2004) combined the WMAP 
data with the three-dimensional power spectrum from over 200,000 galaxies in 
the Sloan Digital Sky Survey (SDSS). Their analysis tightened many of the 
WMAP error bars and gave somewhat lower r e = 0.121qqq. There is still dis- 
agreement between the redshift z r ~ 6 found from (H I) Lya opacity and the 
"early reionization" estimates, z r = 10-20, inferred from the CMB. The optical- 
depth estimates from WMAP are consistent with a wide range in r e , reflecting 
the broad likelihood function (see Figure 4 of Spergel et al. 2003). The analysis 
of second-year WMAP data may yield a more precise determination of r e . 

To elucidate the dependence of r e on the epoch of reionization, we now 
provide the following simplified derivation. Assuming instantaneous, complete 
ionization at redshift z r , we can write r e as the integral of n e aTd£, the electron 
density times the Thomson cross section along proper length. For a standard 
ACDM cosmology, we have (dtjdz) = c{dt/dz) = (1 + z)~^\cj H(z)], where 
H(z) = H Q [n m (l + zf + O a ] 1/2 and Vt m + A = 1. We adopt n e = n H (l + y) 
and nn = [(1 — Y)p CI /mH](l + zf, where Y = 0.244 (He mass fraction) and 
y ~ 0.0807 (He fraction by number). We further assume that helium is singly 
ionized (He II at z > 3). The above integral can be done exactly: 



T P . = 



_c_\ / 2Q b 
HqJ \3Q n 



p cr (l-Y)(l+y)a T 



{n m (l + z r f + Qa} 1 / 2 - 1 



« (0.0376/i) (J^j [{n m (i + z r f + n A } 1/2 - 1] . (2) 

In the above numerical expression, we can substitute the WMAP-concordance 
values (Spergel et al. 2003; Bennett et al. 2003), Q b h 2 = 0.0224±0.0009, 9. m h 2 = 
0.135±g;gg|, and h = 0.7ll^|, where h = [H /(W0 km s" 1 Mpc" 1 ] is the scaled 
Hubble constant. The critical density is p CT = (1.879 x 10~ 29 g cm~ 3 )/i 2 . For 
large redshifts, fl m (l + zf 2> Qa, and the integral simplifies to 



2Q h 



p CT (l-Y)(l + y)a T 



/TM . (0 , ^{l+zrf 2 « (0.00229)(l+z r ) 3 / 2 , 

(3) 

nearly independent of h. Additional contributions to r e arise from residual elec- 
trons (10 -3 - 3 ionized fraction) frozen out after recombination (Seager, Sasselov, 
& Scott 2000) and from partial reionization prior to z r by the first stars (Venkate- 
san, Tumlinson, & Shull 2003) and early X-ray sources (Venkatesan, Giroux, h 
Shull 2001; Ricotti & Ostriker 2004). These additional ionization sources could 
contribute extra electron scattering in the amount Ar e ~ 0.02. We now invert 
the above relation to estimate the reionization redshift, 

/ T P \ 2 / 3 

(1 + *)«(16.2)UM . (4) 



Here, we scaled to r e = 0.15, equal to the WMAP-measured value of optical 
depth, r e = 0.17 ± 0.04, reduced by Ar e = 0.02 arising from additional electron 
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scattering prior to z r . A re-analysis of WMAP data (Tegmark et al. 2004), 
including prior knowledge of the SDSS galaxy power-spectrum, found slightly 
lower optical depth, r e = 0.1241q;qP an d reionization at z r = 14.4jl|;|. The 
large (la) likelihood range allows an optical depth near the value, r e « 0.05, 
expected from from eq. (4), assuming a reionization epoch z r 6 — 7 consistent 
with H I Lya opacity. The agreement would be particularly good if one allows 
for sources of partial ionization at z > z r . 



3. FUSE Observations of He II Absorbers 

Over the last decade, IGM studies have been enhanced by high-resolution optical 
and UV spectra of quasar absorption lines. Studies in H I and He II Lya are 
particularly important for understanding the transition of the high-redshift IGM 
from a neutral to ionized medium. The He II absorption becomes detectable 
in the far ultraviolet at redshifts z > 2.03 with FUSE and at z > 2.79 with 
the Hubble Space Telescope (HST). Moderate-resolution FUSE spectra of He II 
absorption (Kriss et al. 2001; Shull et al. 2004; Zheng et al. 2004) confirm the 
general theoretical picture (Cen et al. 1994; Dave et al. 2001) of the IGM as 
a fluctuating distribution of baryons, organized by the gravitational forces of 
dark matter and photoionized by high-redshift quasars and starburst galaxies 
(Haardt & Madau 1996; Fardal, Giroux, & Shull 1998, hereafter FGS). 

The He II reionization epoch appears to occur at z « 2.8 ± 0.2 (Kriss et al. 
2001; Smette et al. 2002; Shull et al. 2004), based on several sightlines studied 
with HST and FUSE. This effect is demonstrated in Figure 1, which shows the 
patchy regions of broad-band (He II Lya) opacity between z = 2.7 — 2.9 toward 
the QSO HE 2347-4342. Helium reionization is delayed by 1.4 Gyr (Az « 3) 
from hydrogen reionization, probably because He II is more difficult to ionize, 
requiring a harder (4 ryd continuum), probably from AGN. It is also possible that 
partial He II reionization could occur at z > 6, from a population of ultra-hot, 
zero-met allicity massive stars (Tumlinson & Shull 2000; Venkatesan, Tumlinson, 
& Shull 2003). 

FUSE observations of He II can be combined with the corresponding H I 
Lya lines at optical wavelengths, to yield the ratio, rj = N(He II)/N(H I). As we 
explain, measurements of rj probe the shape of the high-redshift metagalactic 
radiation field, presumably from quasars and massive stars. Figure 2 shows a 
portion of the combined observations of He II and H I absorption for the Lya 
forest at z ~ 2.3 — 2.6. The He II absorption is generally much stronger than 
H I by a factor rj =N(He II)/N(H I), predicted theoretically (FGS) to be 50- 
100 for photoionization by quasar backgrounds. For optically thin lines, we can 
write rj ~ 4rH e ii/7"Hi- The greater strength of He II arises because it is harder 
to photoionize than H I, owing to lower fluxes and cross sections at its ionizing 
threshold (hvT = 54.4 eV) and because He III recombines 5 times faster than 
H II. Variations in rj can be used to diagnose the sources and fluctuations of the 
metagalactic ionizing spectrum in the continua of H I (1 ryd) and He II (4 ryd). 
He II absorption is also a good diagnostic of absorption from low-density regions 
("voids") in the IGM. Because THeii = (?7/4)thi (rj S> 1), He II can be tracked 
into much lower density regions than H I. He II absorption can also be used to 
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Figure 1. The epoch of gradual He II reionization at ZHeii ~ 2.8±0.2 toward 
HE 2347-4342 (Shull et al. 2004) is illustrated by the continuum flux transmis- 
sion through He II absorption in several well-defined redshift windows. Note 
the three regions of high transmission between z = 2.7 — 2.8, flanked by low 
transmission regions. The transmission slowly recovers at z < 2.6. 



probe collisionally ionized gas at T rs 10 5 K, produced by shock-heating during 
structure formation (Cen & Ostriker 1999; Dave et al. 2001). 

Following the initial FUSE discovery and analysis paper (Kriss et al. 2001), 
Shull et al. (2004) and Zheng et al. (2004) performed in-depth analyses of the 
He II and H I absorption toward HE 2347-4342, an AGN at redshift z cm = 2.89. 
We probed the IGM at redshifts z = 2.3-2.9, using spectra from FUSE and 
the Ultraviolet- Visual Echelle Spectrograph (UVES) on the VLT telescope. We 
focussed on two major topics: (1) small-scale variability (Az « 10~ 3 ) in the ratio 
r\ = N(He II)/N(H I); and (2) an observed correlation of high-?] absorbers (soft 
radiation fields) with voids in the (H I) Lya distribution. These rj variations 
(Figure 3) probably reflect fluctuations in the ionizing sources on scales of 1 
Mpc, modified by radiative transfer through a filamentary IGM whose opacity 
controls the penetration of 1-5 ryd radiation over 30-40 Mpc distances. In 
photoionization equilibrium, the He II/H I ratio can be expressed: 



_ nHemaHgH r HI _ n -„, • / '" i:i • 

nmi a \n r Heii «/Heii(3 + ai) 



Here, , a^jj, Thi, and THeii are the case-A recombination rate coefficients 
and photoionization rates for H I and He II, and Jhi and JHeii are the specific 
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Figure 2. Overlay of the Lya absorption of He II (top panel - FUSE) and 
of H I (bottom panel VLT/UVES) toward HE 2347-4342. In the VLT data, 
"filament" points are defined as those with thi > 0.05. Wavelengths of H I 
Lya were divided by 4 for alignment. Note the strong He II absorption (top). 



intensities of the radiation field at the H I (912 A) and He II (228 A) edges. The 
parameters a\ and «4 are the local spectral indices of the ionizing background 
at 1 and 4 ryd, respectively, which provide minor corrections to the photoioniza- 
tion rates. We adopt case-A hydrogenic recombination rates to H I and He II, 
appropriate for the Lya forest absorbers with very low neutral fractions. Over 
the temperature range 16,000 K < T < 32,000 K, we approximate the case-A 
recombination rates coefficients as aHeii = (1-36 x 10~ 12 cm 3 s -1 )?^' 700 and 
a m = (2.51 x 10~ 13 cm 3 s" 1 )^ - 755 , where T 4 . 3 = (T/10 4 - 3 K). We assume that 
nHe/^H = 0.08 (Y = 0.244) and that H and He are fully ionized (n e = 1.16rtff) 
after reionization (z < 3). 

Owing to photon statistics and backgrounds, we can measure optical depths 
over the ranges 0.1 < th c ii < 2.3 and 0.02 < thi < 3.9, and we can reliably 
determine values of r) ~ 4rH e ii/THi over the range 0.1 to 460. Values i] = 
20 — 200 are consistent with models of photoionization by quasars with observed 
spectral indices ranging from a s = — 3 (Figure 4). Values rj > 200 may 
require additional contributions from starburst galaxies, heavily filtered quasar 
radiation, or density variations. Regions with r/ < 30 may indicate the presence 
of local hard sources. We find that r\ is higher in "void" regions, where H I is 
weak or undetected and ~80% of the path length has rj > 100. These voids 
may be ionized by local soft sources (dwarf starbursts) or by QSO radiation 
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Figure 3. Observed wide variations in the ratio r\ = N(He II)/N(H I) in 
the IGM toward HE 2347-4342 (Shull ct al. 2004). These variations imply 
substantial fine-grained fluctuations in the ionizing radiation field at 1 ryd 
(H I) and 4 ryd (He II) on scales of Az = 0.001 or 1 Mpc. 



softened by escape from the AGN cores or transfer through the "cosmic web". 
The apparent differences in ionizing spectra may help to explain the 1.4 Gyr lag 
between the reionization epochs of H I (zm ~ 6.1 ± 0.3; Gnedin 2004) and He II 
(zHeii ~ 2.8 ± 0.2; Kriss et al. 2001; Shull et al. 2004). 

Many of the variations in the ionizing radiation field appear to come from 
intrinsic differences in QSO spectra (Fig. 4), found by direct measurements of 
the redshifted EUV continua using HST (Telfer et al. 2002) and FUSE (Scott 
et al. 2004). Figure 5 shows some of the actual AGN spectra on which these 
conclusions are based. It is apparent that even low-redshift QSOs and Seyferts 
have wide differences in their EUV continuum slopes, which may be luminosity 
dependent (Scott et al. 2004). In addition to intrinsic source variations, the 
metagalactic radiation field is heavily filtered and reprocessed by intervening 
IGM (see FGS). Because the He II opacity is so strong (77 S> 1), fluctuations at 
4 ryd are probably larger than those at 1 ryd. Unfortunately, we must conclude 
that the sources, spectra, and fine-grained spatial variations in the metagalactic 
ionizing radiation field are still not well understood. 

Future studies of He II should proceed along several fronts. First and fore- 
most, we need to find additional high-redshift AGN targets for which we can 
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Figure 4. Observed distribution of the observed QSO ionizing spectral in- 
dices, a E uv, for radio-quiet AGN observed by HST (Telfer et al. 2002). Other 
curves show the spectral indices a s , with J„ oc v~ as , that are required to 
reproduce the observed values of He II/H I ratios, using eq. [5]. Note the sys- 
tematically softer radiation fields (high 77, large a s ) in the voids (thi < 0.05) 
compared to filaments (thi > 0.05). The breadth of these distributions and 
the clear offset between voids and filaments suggest that additional spectral 
softening may be present. 



repeat the He II/H I measurements and rj analysis. Quasar H1700+6416 at 
z cm = 2.74 is one such FUSE target, but more are needed to identify the cosmic 
variance in z r (He II). Additional theoretical work is needed, using cosmological 
simulations that include radiative transfer of the 1 ryd and 4 ryd radiation in- 
tensities and their spatial fluctuations. Radiation in the 3-4 ryd band may be 
especially important for the photoionization corrections to C III, C IV, Si III, 
and Si IV on which IGM metallicities are based. Finally, it would be helpful 
to characterize the mean intrinsic EUV spectra of AGN, at z < 0.5, which are 
still highly uncertain (Figure 5) despite heroic efforts with both HST and FUSE. 
There is currently only a small wavelength window at which the rest-frame EUV 
is observable, and it does not extend to A < 228 A, the 4 ryd continuum that 
controls He II ionization. 
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Figure 5. Normalized distribution of fits to the rest-frame UV and EUV 
continuum spectra of 73 AGN (Clausen & Shull 2004), taken from FUSE data 
analyzed and fitted by Scott et al. (2004). We have shifted these data to the 
rest-frame, K-corrected the UV fluxes, and normalized them to unity at 1 ryd 
(912 A). The solid lines show actual data, and dotted lines are extrapolations. 
Note the wide variations in spectral indices at A = 650 — 900 A. 
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